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ABSTRACT 

We present Hubble Space Telescope imaging and spectroscopic observations 
of three Brightest Cluster Galaxies, Abell 1836-BCG, Abell 2052-BCG, and 
Abell 3565-BCG, obtained with the Wide Field and Planetary Camera 2, the 
Advanced Camera for Surveys and the Space Telescope Imaging Spectrograph. 
The data provide detailed information on the structure and mass profile of the 
stellar component, the dust optical depth, and the spatial distribution and kine- 
matics of the ionized gas within the innermost region of each galaxy. Dynamical 
models, which account for the observed stellar mass profile and include the contri- 
bution of a central supermassive black hole (SBH), are constructed to reproduce 
the kinematics derived from the Ha and [N 11]AA6548,6583 emission fines. Se- 
cure SBH detection with M. = 3.611°;^^ x 10^ Mq and M. = IMt^'^ x 10^ 
M0, respectively, are obtained for Abell 1836-BCG and Abell 3565-BCG, which 
show regular rotation curves and strong central velocity gradients. In the case 
of Abell 2052-BCG, the lack of an orderly rotational motion prevents a secure 
determination, although an upper limit of M, < 4.60 x 10^ Mq can be placed on 
the mass of the central SBH. These measurements represent an important step 
forward in the characterization of the high- mass end of the SBH mass function. 

Subject headings: black hole physics — galaxies: elliptical and lenticular, cD — 
galaxies: nuclei — galaxies: kinematics and dynamics 



1. Introduction 

Within the past decade, the focus in supermassive black holes (SBHs) studies has moved 
from the dynamical measurement of SBH masses, M,, in nearby galaxies (see review by Fer- 
rarese & Ford 2005), to the characterization of scaling relations connecting M, to the large 



^ Based on observations made with ESO telescopes at the La Silla Paranal Observatory under programme 
ID 279.B-5004(A). 
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scale properties of their hosts (Kormendy & Richstone 1995; Ferrarese & Merritt 2000; Geb- 
hardt et al. 2000; Graham et al. 2001; Ferrarese 2002; Marconi & Hunt 2003). Such 
relations, combined with the knowledge of the galaxy luminosity or velocity dispersion func- 
tions, lead to a direct determination of the local SBH mass function and, by comparison with 
the energetics of high redshift AGNs, accretion history (e.g. Marconi et al. 2004; Shankar 
et al. 2004; Benson et al. 2007; Tundo et al. 2007; Graham et al. 2007; Lauer et al. 2007). 
Furthermore, the tightness of the relations linking galactic properties to M, is indicative of a 
formation/evolutionary history in which SBHs and galaxies are causally connected. Indeed, 
feedback from AGN activity is believed to play an important, perhaps even fundamental 
role in the evolution of galaxies (e.g. Binney & Tabor 1995; Suchkov et al. 1996; Ciotti & 
Ostriker 2001; Schawinski et al.2006; Springel et al. 2005; Hopkins et al. 2007; Di Matteo 
et al. 2007). 

In this context. Brightest Cluster Galaxies (BCGs) are of particular interest. Their 
privileged location at the center of a massive cluster implies that they have undergone a 
particularly extensive merging history (Khochfar & Silk 2006) and are likely to host the 
most massive SBHs in the local Universe (Yoo et al. 2007). The latter conclusion is also 
supported by all scaling relations which are known to be obeyed by local SBHs, that predict 
the most massive SBHs to reside in most massive galaxies. As such, BCGs constitute an 
excellent laboratory to search for the local rehcs of the most powerful high-redshift quasars 
(Willott et al. 2003; Vestergaard 2004), and to investigate the role of mergers in the black 
hole mass function. 

Unfortunately, direct dynamical measurements of SBHs masses in BCGs are exceedingly 
difficult - with only two such measurements made to-date (in M87 and NGC 1399, Harms et 
al. 1994; Macchetto et al. 1997; Houghton et al. 2006). The reason for this is simple: SBH 
mass measurements based on resolved kinematic tracers (gas or stars) need to be carried 
out within the SBH "sphere of influence", i.e. the region of space within which the SBH 
dominates the overall gravitational potential. For a 10^ Mq SBH, the sphere of influence 
becomes unresolved for optical measurements with the Hubble Space Telescope at distances 
beyond ~ 100 Mpc (sec Figure 44 of Ferrarese & Ford 2005). Few BCGs are located within 
this limit. To compound the problem, bright ellipticals are characterized by shallow stellar 
density profiles and faint central surface brightnesses (e.g., Ferrarese et al. 1994; Lauer et 
al. 1995; Rest et al. 2001; Ferrarese et al. 2006) making the detection of stellar absorption 
lines with HST prohibitively expensive. The latter difficulty can be overcome with the use 
of gas dynamics. Emission lines (mainly Ha and [N II]) from the gas are bright and easily 
measured. If the gas is confined in a disk, there is little ambiguity in the velocity distribution 
(Ho et al. 2002), and since the method was first used (Harms et al. 1994; Ferrarese et al. 
1996; Ferrarese & Ford 1999) an increasing amount of attention has been devoted to the 
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theoretical aspects of the dynamical modeling (Maciejewski & Binney 2001; Barth et al. 
2001; Cappellari et al. 2002; Coccato et al. 2006). 

The lack of a secure characterization of the SBH mass function above the 10^ Mq 
mark is troublesome. Lauer et al. (2007) suggest that the relations between M, and host 
bulge luminosity, Lbuige, (Kormendy & Richstone 1995) and velocity dispersion (Ferrarese & 
Merritt 2000; Gebhardt et al. 2000) would predict significantly different M, if extrapolated 
above 10^ Mq. In particular, the M, — a relation would predict less massive SBHs in BCGs 
than the M, — Lbuige relation, due to the slower increase of a with galaxy luminosity observed 
for BCGs compared to the bulk of the Early- Type galaxy population (Bernardi et al. 2007). 
The difference in the predicted mass is such to significantly affect (by an order of magnitude) 
the high-mass end of the local SBH mass function, von der Linden et al. (2007) argue that 
the shallower dependence of a on L applies to BCG but not to comparably massive non-BCG 
galaxies, which instead follow the canonical Faber- Jackson relation defined by less massive 
systems. This implies that BCGs and non-BCG galaxies of comparable mass must occupy 
a different locus in either, or both, the M, — a and M, — Lbuige- This result is in contrast 
with the findings of Batcheldor et al. (2007), who argue that SBHs masses predicted from 
NIR luminosities (Marconi & Hunt 2003) are consistent with those predicted from a. They 
attributed the discrepancy noted by Lauer et al. (2007) (who used l^— band luminosities) to 
a bias introduced by the inclusion, when computing Lbuige, of extended blue envelopes around 
BCGs. At present these ambiguities prevent us from testing any theoretical prediction on 
the high-mass end of the black hole mass function based on the observed quasar abundances 
and merger histories. 

In this paper, we analyze the kinematics of the ionized gas in the nuclear region of three 
BCGs - Abell 1836-BCG, Abell 2052-BCG, and Abell 3565-BCG - in order to constrain the 
mass of the central SBHs. The data were obtained using the Space Telescope Imaging Spec- 
trograph (STIS) on board the Hubble Space Telescope (HST), supplemented with imaging 
data from the Advanced Camera for Surveys (ACS) and the Wide Field and Planetary Cam- 
era 2 (WFPC2). Furthermore, the central stellar velocity dispersion of Abell 1836-BCG was 
obtained from ground based spectroscopy. The paper is organized as follows. The criteria of 
galaxy selection are presented in §21 ACS and STIS observations are described and analyzed 
in §[|S]and|ll respectively. Ground based spectroscopic observations of Abell 1836-BCG are 
described and analyzed in ^ 

The SBH masses of Abell 1836-BCG and Abell 3565-BCG, and an upper limit for the 
SBH mass of Abell 2052-BCG are derived in ^ In §7]results are compared to the predictions 
of the SBHs scaling laws and conclusions are given. 
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2. Galaxy sample 

Naturally, the place to search for the most massive black holes in the universe is in 
the most nearby and most massive BCGs, where the problems of angular resolution to 
measure black hole masses are minimized. The galaxies discussed in this contribution were 
selected from the BCG sample observed by Laine et al. (2003) with the HST Wide Field and 
Planetary Camera 2 (WFPC2). Dynamical studies in giant ellipticals are performed more 
effectively by modeling gas, rather than stellar kinematics (see discussion in Ferrarese & Ford 
2005), both because of the low stellar surface brightness that characterizes the core regions 
of such systems (e.g. Ferrarese et al. 2006) and because when applied to slowly-rotating, 
pressure supported galaxies, stellar dynamical models depend heavily on the knowledge of the 
anisotropy in the stellar velocity dispersion, a difficult quantity to constrain observationally 
(e.g. van der Marel 1998). Eleven of the 81 galaxies in the Laine et al. (2003) sample appear 
to contain well defined nuclear dust structures, which can trace regular gas kinematics (Ho et 
al. 2002). After selecting the most regular nuclear dust structures among those eleven objects 
and after excluding galaxies for which existing ground based spectra showed no evidence of 
nuclear Ha and [N II] AA6548, 6583 in emission, we selected the four BCGs predicted to host 
the most massive SBHs (based on the correlations between black hole mass, M,, and the host 
bulge stellar velocity dispersion, a, and/or S-band luminosity. Lb) and for which the SBH 
sphere of infiuence could be well resolved using the 0".l wide slit of HST's Space Telescope 
Imaging Spectrograph (STIS). 

One galaxy (Abell 2593-BCG) was not observed due to the premature failure of STIS. 
The sample discussed in this paper therefore includes three BCGs, namely Abell 1836-BCG 
(PGC 49940), Abell 2052-BCG (UGC 9799) and Abell 3565-BCG (IC 4296). The main 
properties of these galaxies are summarized in Table [H 

3. Observations, data reduction, and analysis: Imaging 

3.1. Observations and data reduction 

The luminosity density and ionized gas distribution of the target galaxies, both of which 
are required for the dynamical modeling described in ^ have been constrained using HST 
ACS and WFPC2 images. Basic information about the data can be found in Table [2J 

The ACS images were obtained with the High Resolution Channel (HRC) as part of 
program GO-9838 (P.l. L. Ferrarese). The HRC detector is a 1024 x 1024 SITe CCD with a 
plate scale of 0'f028 x 0'f025 pixel"^. Because of the large, but well characterized, geometric 



- 6 - 



distortion affecting the ACS, the 29" x 26" HRC field of view projects on the plane of the 
sky as a rhomboid with x and y axes forming a 84°2 angle. Each galaxy was observed with 
two filters: F435W (which approximates the Johnson 5-band) and the narrow-band ramp 
filter FR656N, tuned to cover a 130 A wide spectral region centered on the redshifted Ha 
and [N II] AA6548, 6583 emission expected at the nuclear location. To help in identifying and 
correcting cosmic ray events, two back-to-back exposures were taken with each filter. 

The WFPC2 images, all employing the F814W filter (similar to the Kron-Cousin / 
band), belong to programs GO-5910 (P.I. T. Lauer) and GO-8683 (P.I. R.P. van der Marel) 
and were downloaded from the public HST archive. Two images with a total exposure time 
of 1000 s are available for each Abell 1836-BCG and Abell 2052-BCG, while ten images are 
available for Abell 3565-BCG, for a total 11,600 s exposure. In all cases, the galaxies are 
centered on the higher resolution Planetary Camera (PC), which consists of an 800 x 800 
pixels Loral CCD with a plate scale of 0"0455 pixel"^, yielding a nominal field of view of 
36" X 36". For both ACS and WFPC2 images, the telescope was guiding in fine lock, giving 
a typical rms tracking error per exposure of 0"005. 

All images were calibrated using the standard reduction pipelines maintained by the 
Space Telescope Science Institute (PyRAF/CalACS for the ACS images and IRAF/CalWFPC 
for the WFPC2 images). Reduction steps include bias subtraction, dark current subtrac- 
tion and fiat-fielding, as described in detail in the ACS and WFPC2 instrument and data 
handbooks (Raggett et al. 2002; Heyer et al.2004; Pavlovsky et al. 2004a,b). 

Subsequent analysis was performed using standard tasks in For each galaxy, the 

alignment of images obtained with the same instrument and filter was checked by comparing 
the centroids of stars in the field of view; in all cases but one, the alignment was found to 
be better than a few hundredths of a pixel and the images were combined without further 
processing using IRAF/IMCOMBINE. The exception are the WFPC2 images of Abell 3565- 
BCG, which were dithered by ~ 1" relative to one another. These images were aligned using 
IMSHIFT prior to being combined. 

In combining the images, pixels deviating by more than three times the local sigma - 
calculated from the combined effect of Poisson and read-out noise - were flagged as cosmic 
rays and rejected. Finally, the resulting ACS images were corrected for geometric distortion 
using PYRAF/PYDRIZZLE. 



^IRAF is distributed by the National Optical Astronomy Observatories which are operated by the Asso- 
ciation of Universities for Research in Astronomy (AURA) under cooperative agreement with the National 
Science Foundation. 
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The final, calibrated ACS and WFPC2 images of each galaxy are shown in Figures [T]l3] 
andlH respectively. 



Isophotal parameters (coordinates X, Y of the isophotal center; surface brightness profile 
/((/)); ellipticity e; major axis position angle 6; and deviations of the isophotes from pure 
ellipses) were measured using the IRAF task ELLIPSE (Jedrzejewski 1987). For each semi- 
major axis length, the best fitting set of parameters X, Y, e and 6 are those that minimize 
the sum of the squares of the residuals between the data and the first two moments of the 
Fourier series expansion of the azimuthally-sampled intensity 



The third and fourth order moments (^3, B^, Ai^ and B^^ describe three and four- 
fold deviations of the isophotes from pure ellipses, respectively. The -B4 term in particular 
describes the shape by distinguishing boxy [B/^ < 0) from disky {B^ > 0) galaxies. 

The isophotal semi-major axis was increased logarithmically, with each isophote being 
fitted at a semi-major axis 10% longer than that of the isophote preceding it. The dust 
structure in Abell 2052-BCG is localized enough not to affect significantly the recovery of 
the underlying stellar distribution, so all parameters were allowed to vary freely while fitting 
the isophotes. For the other galaxies, a fully unconstrained isophotal solution could be found 
only in the regions unaffected by dust (r > 0".6 and r > 0''9 for Abell 1836-BCG and Abell 
3565-BCG, respectively). In the areas affected by dust, the isophotes were fitted by fixing 
both ellipticity and position angle to the average values obtained in the dust-unaffected 
region ( 0".6 < r < 10" and O'^Q < r < 10" for Abell 1836-BCG and Abell 3565-BCG, 
respectively). In all cases, pixels deviating by more than three times the standard deviation 
of the mean intensity along each trial isophote were excluded from the fit, thereby avoiding 
contamination by foreground stars, companion galaxies (as in the case of Abell 2052-BCG), 
globular clusters, and bad pixels. All isophotal parameters are plotted in Figures [5H7] as a 
function of the "geometric mean" radius r^eo, defined as a[l — e{a)]^^'^, a being the semi-major 
axis length. 

Conversion of /(0) (which is given in counts per seconds) to surface brightness in the 
AB magnitudes system was performed following Holtzman et al. (1995) for the WFPC2 
images, and Sirianni et al. (2005) for the ACS images: 



3.2. Isophotal analysis 




(1) 



k 
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/^F435W 



■2.5 log/,((/))F8i4W + 14.998, 
2.5 logJ,(0)F435W + 16.150, 



(2) 
(3) 



where /s(0)f814W and /s(0)f435w represent the azimuthally averaged, background sub- 
tracted intensity (in counts per second per pixel) in F814W and F435W respectively. Because 
the galaxies completely fill the field of view of the ACS/HRC and WFPC2/PC detectors, a 
nominal sky background, as listed in the ACS (Pavlosky et al. 2004a) or WFPC2 (Heyer et 
al. 2004) Instrument Handbooks, was adopted. 



Correction for dust obscuration is necessary to recover the intrinsic stellar luminosity 
density and the spatial distribution of the ionized gas within each galaxy. We adopted the 
procedure described in Ferrarese et al. (2006), which works under the assumption that 
the F435W-F814W intrinsic color of the galaxy in the regions obscured by dust can be 
estimated by linearly interpolating across the dust area (or extrapolating, if the dust affects 
the center) the azimuthally- averaged color measured in the regions unaffected by dust. The 
dust absorption at each pixel then follows by comparing the intrinsic and observed (extincted) 
color maps, once the ratio of the absorption in F435W and F814W is known (see equation 
2 of Ferrarese et al. 2006). The latter, Af435wMf814W = 2.2814 (where Ax is given in 
magnitudes), has been derived using the extinction law of Cardelli et al. (1989), integrated 
over the filter passbands, and assuming a ratio of total-to-selective absorption Ry = 3.1 (see 
also Schlegel et al. 1998). It is further assumed that the dust lies in the galaxy foreground 
(the "screen" approximation); the values quoted below for the dust absorption therefore 
represent firm lower limits to the true absorption. Once the absorption in a specific passband 
is known, the intrinsic fiuxes are easily recovered simply by correcting the measured fiuxes 
for the derived magnitude loss due to dust obscuration. 

In producing the F435W-F814W color maps, the ACS images were rotated (using the 
information about the position angle of the telescope axis contained in the header) and 
resampled using PYDRIZZLE in order to match the orientation and scale of the lower-resolution 
WFPC2 images. 

Dust absorption maps are shown in the left panel of Figures [5HTU1 In Abell 1836- 
BCG, the average extinction (corrected to the banclfl) within the regular nuclear dust 
disk is {A{V)) = 0.27 mag, with values as high as A(y)max = 0.85 mag (in the screen 
approximation). The nuclear dust disk in Abell 3565-BCG is quite similar, with A{V)max = 



3.3. Dust obscuration and optical depth map 
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1.1 mag and {A{V)) = 0.21 mag. Recovering a dust extinction map for Abell 2052-BCG, in 
which the dust is distributed in an irregular filamentary pattern and the amount of extinction 
is more modest, was more challenging. Within the inner ^ 1", we find A{V)max = 0.35 mag 
and {A{V)) = 0.14 mag. At larger radii, dust patches with lower levels of dust obscuration, 
although visible in the color maps, could not be properly modeled. These areas, which affect 
less than 50% of pixels along any given isophote, were simply masked when deriving the 
galaxy's surface brightness profile. 

The extinction corrected F814W surface brightness profiles of the sample galaxies are 
shown in Figures [5H71 



The deprojected luminosity density, T{x,y,z), needed to constrain the stellar gravita- 
tional potential, was recovered from the redder (and thereby less affected by dust obscu- 
ration), extinction corrected, WFPC2/F814W images, using the multi-Gaussian expansion 
(MGE) method (Monnet at al. 1992; Emsellem et al. 1994), as implemented by Cappellari 



Briefly, the MGE procedure starts by determining the galaxy's average ellipticity, po- 
sition angle, and coordinates of the luminosity-weighted center. The galaxy image is then 
divided in four quadrants, and a number of photometric profiles is measured along sectors 
uniformly spaced in angle from the major axis to the minor axis. Surface brightness profiles 
from sectors in the four quadrants are averaged together (an acceptable simplification for 
our galaxies since position angle and isophotal center do not show a significant radial depen- 
dency), and each is then fitted as the sum of Gaussian components. The best-fitting MGE 
model surface brightness is then determined iteratively by comparison with the observed 
surface brightness, after having been convolved with the instrumental PSF (generated, in 
our case, using the TinyTIM package, Krist & Hook 1999, and parametrized as the sum of 
circular Gaussians). The Gaussian width coefficients of the MGE model were constrained to 
be a set of logarithmically-spaced values, thus simplifying the fitting algorithm into a general 
non-negative least-squares problem for the corresponding Gaussian amplitudes. During the 
entire procedure, pixels for which correction for dust obscuration did not perform correctly 



■^Extinction in all relevant passbands are related to A(V) as v4f435w/^v — 1-362, Apsuw/^v — 0.597, 
^fr656nMv = 0.781,0.783, and 0.805 for Abell 1836-BCG, Abell 2052-BCG, and Abell 3565-BCG, respec- 
tively (Cardelli et al. 1989; Schlegel et al. 1998). 

''The code is publicly available at [http: / /www. strw.leidenuniv.nl/~nicappell /idl / 



3.4. Deprojected stellar density models 
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(as in the case of Abell 2052-BCG) were masked. 

The best-fitting MGE model (prior to PSF convolution) is then deprojected to a lu- 
minosity density r(r) under the assumption that the galaxy is spherically symmetric - a 
justifiable simplification since the isophotes are nearly circular throughout the main body of 
our galaxies, including the radial range where the ionized-gas kinematics probes the galaxy 
potential (where however uncertainties are larger due to the presence of dust residuals, see 
Figures [3H7]). For a radially-independent mass-to-light ratio, (M/L)^, the stellar mass density 
p(r) is then simply expressed as p(r) = (M/L)j,r(r), and the stellar gravitational potential 
can be expressed in terms of error functions. The best-fitting, PSF-convolved MGE models 
to the observed surface-brightness profiles are shown in Figures [TTHT3I along with the depro- 
jected luminosity density profiles (prior to PSF-convolution), and the corresponding circular 
velocity curves [for {M/L)^= 1]. 



3.5. Ionized-gas distribution 

Continuum-free emission-line images for the sample galaxies were obtained by subtract- 
ing the WFPC2/F814W images from the ACS/FR656N images which isolate the spectral 
region characterized by the redshifted Ha and [N II] AA6548, 6583 emission lines. All im- 
ages were previously extinction corrected as described in § 13. 3^ after having been resam- 
pled and rotated in order to match the pixel scale and orientation of the lower-resolution 
WFPC2/F814W images. Before subtraction, the FR656N images were multiplied by a fac- 
tor equal to the mean ratio between the F814W and resampled FR656N images, calculated 
in emission-free regions. The resulting emission line maps are shown in the central panels 
of Figures [SHTUl The intrinsic (prior to PSF broadening) surface-brightness distribution of 
the ionized gas, a key input to our modeling, was obtained independently using both the 
iterative method based on the Lucy-Richardson algorithm (Richardson 1972; Lucy 1974) 
and the MGE method already described in § 13. 4[ According to van den Bosch & Emsellem 
(1998) the two methods give consistent results, however, one might be preferable to the 
other depending on the specifics of the images (e.g., signal-to- noise ratio S/N, presence of 
asymmetric features, etc.). Both methods employed the PSF generated using TINY TIM 
for the WFPC2/F814W images. The Lucy- Richardson algorithm, as implemented in the 
IRAF task LUCY, was found to give satisfactory results in the case of Abell 2052-BCG and 
Abell 3565-BCG, for which the method converged after only five and six iterations, respec- 
tively. In the case of Abell 1836-BCG, for which the emission-line images have lower S/N, 
LUCY produced unacceptable noise amplification before convergence could be reached, and 
the intrinsic emission-line distribution was therefore recovered using the MGE method. We 
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note that the latter failed to produce acceptable results for Abell 2052-BCG and Abell 3565- 
BCG, for which the distribution of the ionized gas is less symmetric, violating the four-fold 
symmetry required by the MGE algorithm. The deconvolved surface-brightness distributions 
of the emission line gas are shown in the right panels of Figures [5HTU1 

4. Observations, data reduction and analysis: STIS spectroscopy 

The spectroscopic observations of the sample galaxies were carried out with STIS be- 
tween 2003 August 16 and 2004 February 21 as part of program GO-9838. The G750M 
grating was used at the prime tilt in combination with the 0''lx52" slit for Abell 1836-BCG 
and Abell 2052-BCG, and the 0':2x52" slit for Abell 3565-BCG. The slit width was selected 
to ensure that the sphere of influence of the SBH estimated to reside in each galaxy (based 
on the M, — a and M,— relations) could be comfortably resolved, while at the same 
time allowing for reasonable throughput. The G750M grating covers the wavelength range 
from 6480 A to 7050 A, which includes the Ha, [N II] AA6548, 6583, and [S II] AA6716, 6731 
emission lines. The dispersion is 0.56 A pixel^^, and the instrumental resolution is 0.83 A 
(FWHM, corresponding to (Tinstr ^ 16.0 km s"^ at [N II] A6583) for the 0'flx52" slit, and 
1.5 A ((Tinstr ^ 30.5 km s"^ ) for the 0"2x52" slit. The spatial scale of the 1024x1024 SITe 
CCD is 0'f05 pixel-^ 

For each target galaxy, spectra were taken at three different slit positions, the first 
crossing the nucleus of the galaxy, and two additional ones displaced by one slit width (O'fl 
for Abell 1836-BCG and Abell 2052-BCG, and 0'f2 for Abell 3565-BCG) on either side, 
perpendicularly to the slit axis. All slits were oriented along the photometric major axis. 
In the case of Abell 3565-BCG, for which the wider slit allowed for shorter exposure times, 
one spectrum was taken at each slit position. For Abell 1836-BCG and Abell 2052-BCG two 
spectra were taken at each location, the second shifted along the slit axis by four pixels to 
facilitate removal of bad pixels and other CCD defects. A 90 s tungsten lamp fiat exposure 
was taken after each galaxy spectrum, to correct for the fringing affecting the STIS CCD at 
wavelengths longer than 7000 A. A log of the STIS observations can be found in Table [31 

The spectra were reduced using the IRAF/CALSTIS reduction pipeline maintained by 
STScI. The basic reduction steps included overscan subtraction, bias subtraction, dark sub- 
traction, flatfield correction, wavelength calibration as well as correction for geometric dis- 
tortion. Correction for fringing was performed with the IRAF/DEFRINGE task. In the case 
of Abell3565-BCG, for which only one spectrum is available at each slit location, cosmic 
ray events and hot pixels were removed using the task LACOS_SPEC (van Dokkum 2001). 
For the other two galaxies, the two spectra obtained at the same slit position were aligned 
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with IRAF/IMSHIFT using the center of the stellar-continuum radial profile as a reference, 
and then combined using IRAF/IMCOMBINE. Cosmic ray rejection was performed with the 
CRREJECT option in IMCOMBINE, additionally, pixels fiagged by running LACOS_SPEC sepa- 
rately on each spectrum were also removed. The sky background was determined from the 
combined, calibrated spectra between 15" and 30" from the nucleus and then subtracted. 
Finally, individual spectra were extracted every 0"05 (corresponding to 1 pixel) up to ~ 0"2, 
0':3, and 0"7 from the center for Abell 1836-BCG, Abell 2052-BCG, and Abell 3565-BCG, 
respectively. At larger radii, the decreasing intensity of the emission lines required binning in 
the spatial direction (see Table [Hand Figures [T^l - [T6l for details) to obtain a S/N of at least 10 
at the peak of the [N II] A6583 emission line in the combined spectrum, a condition deemed 
necessary for an accurate measurement of the line kinematics (§4.1). The final spectra along 
the major and offset axes of the sample galaxies are plotted in Figures fT^b)-[T6l(b) . 

The location of the slit relative to the galaxy center was checked using STIS images taken 
after acquisition of the target galaxies. A series of "synthetic" slits were extracted from the 
acquisition images by averaging four adjacent columns for Abell 1836-BCG and Abell 2052- 
BCG (corresponding to a 0"l-wide slit), and eight adjacent columns for Abell 3565-BCG 
(corresponding to a 0"2-wide slit) at the nominal slit orientation. The radial brightness 
profiles in each synthetic slit were then compared to the profiles obtained by collapsing 
the spectrum at the nuclear location along the wavelength direction. The best match was 
determined by minimizing the of the difference between the light profile of the spectrum 
and the light profile extracted from the acquisition image. In all cases, the slits were found 
to be all positioned at the nominal locations. 



4.1. Measurement of the emission lines 

The ionized-gas kinematics were measured from the narrow [N II] A6583 line, the bright- 
est of the lines detected in our spectra. This line is preferred to Ha since the latter is 
more severely blended, and could be significantly affected by emission from circumnuclear 
star- forming regions (e.g., Verdoes Kleijn et al. 2000; Coccato et al. 2006). 

The emission lines' central wavelengths, FWHMs, and intensities were measured fol- 
lowing Beifiori et al. (2008). The best-fitting Gaussian parameters were derived using a 
non-linear least-squares minimization based on the robust Levenberg-Marquardt method 
(e.g.. Press et al. 1992) implemented by More et al. (1980). The actual computation 
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was done usine; the MPFIT algorithirO implemented by C. B. Markwardt under the IDL|j 
environment. 

Because the hne profiles are assumed to be Gaussian, it was often necessary to include 
two components in order to reproduce the extended wings observed for several of the emission 
lines (most notably Ha but, in some cases, also [N II] AA6548, 6583 and [S II] AA6716, 6731). 
Lines in the same doublet ([N II] and [S II]) were assumed to have the same velocity and 
velocity width. A flux ratio of 1:2.967 was assumed for the two [N II] lines, as dictated 
by atomic physics (e.g., Osterbrock 1989). The stellar continuum was approximated with a 
low-order polynomial. 

The major axis spectrum of Abell 1836-BCG displays both a narrow and a broad Ha 
line within the inner 0".25. The broad-line emission appears to be spatially extended, with 
FWHM of 0".17, compared to the typical FWHM of the STIS PSF (~ 0':075). Only the 
—0".l offset spectrum shows a broad Ha component within the inner 0".l (two rows), which 
possibly originates in the same nebular complex that gives rise to the broad-line emission 
observed along the central slit. 

Only the major axis spectrum of Abell 3565-BCG displays broad-line emissions. Within 
the inner 0".15 the spectrum shows an Ha component while from 0"15 to 0"A1 it shows broad 
components in Ha, [N II] AA6548, 6583, and [S II] AA6716, 6731. The broad Ha emission is 
spatially extended, with a FWHM of 0';i3. 

The major-axis and the —0".l offset spectra of Abell 2052-BCG show both narrow and 
broad components in Ha, [N II] AA6548, 6583, and [S II] AA6716, 6731 within the inner 0".30. 
The O'^l offset spectrum shows all the broad components in the inner 0".l. The broad lines 
have a FWHM of O'flO for Ha, 0".12 for the [N II] doublet, 0':i8 for [S II] A6716, and 0'!08 
for [S II] A6731 ([S II] broad lines have the same radial extension). They appear to be less 
extended than for the other two galaxies, although still marginally resolved. 

Figures [T71 - [T9l show the results of the multi-Gaussian fitting for the spectra extracted at 
the nuclear location, and at a location displaced, along the major axis, by ~ 70 parsec (0'fl5, 
O'flO, and 0"25 for Abell 1836-BCG, Abell 2053-BCG, and Abell 3565-BCG respectively). 

Line-of-sight heliocentric velocities, velocity dispersions and line intensities measured 
from the [N II] A6583 line along the major and offset axes of the nuclear dust disks are 
plotted in Figures [T W c)- [W c) and listed in Table HI 



^The updated version of this code is available on http://cow.physics.wisc.edu/~craigm/idl/idl.html 
^Interactive Data Language is distributed by Research System Inc. 
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4.2. lonized-gas kinematics and dust morphology 

In the remainder of this Section, we will comment briefly on the kinematics of each 
individual galaxy. Before doing so, we note that both Abell 1836-BCG and Abell 3565-BCG 
display regular and symmetric velocity flelds, as suggested by the smoothness and regularity 
of the ionized-gas disks. In contrast, Abell 2052-BCG, which is characterized by a more 
irregular emission-line morphology, shows irregular kinematics. 

Abell 1836-BCG. The rotation curve, velocity dispersion proflle, and flux proflle along 
the major axis are relatively regular (Figure [T^). which makes this object appealing for 
dynamical modelling. The rotation velocity increases rapidly within the inner 0''1, reach- 
ing approximately 200 km s~^, and flattening at larger distances. The velocity dispersion 
peaks at about 370 km s~^ at the nuclear location, and declines on either side, although not 
symmetrically. The off-axis kinematics show a similar pattern. 

Abell 2052-BCG. The dust distribution in this galaxy is somewhat irregular, with several 
filaments distributed in a spiral-like pattern (Figure [2]). This lack of symmetry is reflected 
in the flux distribution of the ionized gas, which is also not symmetric with respect to the 
center (lower panels of Figure [T5b). and in the [N II] A6583 kinematics, which could be 
measured out to about 0".3 from the center along each slit positions (Figure At the 

nuclear location, the maximum rotation velocity (200 km s~^) is reached at 0'f05, with the 
velocity decreasing on either side, while the velocity dispersion curve is more symmetric, 
reaching «i 500 km s~^at the nuclear location. The spectra flanking the nuclear spectra are 
characterized by similar velocity dispersion proflles, but very different rotation curves; in 
particular, almost no rotation is observed along the slit displaced by 0"1 to the South-East. 
Because of this lack of regularity in the gas kinematics, we are only able to place an upper 
limit on the SBH mass of Abell 2052-BCG f^ jOj) . 

Abell 3565-BCG. The [N II] A6583 kinematics are measured within 0".8 from the center 
along all three slit locations (Figure [TBT c)). Although the smooth and well deflned dust 
disk in Abell 3565-BCG starts to tilt slightly at large radii from the center (Figure [3]), such 
a warp is clearly visible, and presumably kinematically important, only beyond the region 
where the ionized-gas kinematics is measured. The major-axis rotation velocity has a steep 
gradient, increasing to approximately 250 km s~^ within 0".l. It remains approximately 
constant out to the last observed point along the approaching side, while it decreases to 
about 100 km s~^ for radii larger than 0"5 along the receding side. The velocity dispersion 
profile is fairly symmetric and strongly peaked, reaching 450 km s~^ at the center. The 
offset axes show rather similar kinematics, although the line fluxes are rather different. The 
rotation velocity increases almost linearly with radius, up to ~ 200 km s~^ at ~ 0".3 from 
the center, while the velocity dispersion remains fairly constant at 120 km s~^. Like Abell 
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1836-BCG, Abell 3565-BCG is a promising candidate for dynamical modeling. 

5. Observations, data reduction and analysis: Ground-based spectroscopy 

Although not used in the dynamical analysis, large scale stellar kinematics are necessary 
to place measured SBH masses on the M,—ac relation. Velocity dispersions measurements 
exist in the literature for Abell 2052-BCG and Abell 3565-BCG, but not for Abell 1836-BCG. 
Ground-based spectroscopic observations of Abell 1836-BCG were therefore carried out in 
service mode with the 3.5-m New Technology Telescope (NTT) at the European Southern 
Observatory (ESO) in La Silla (Chile) on 2007 May 20 and July 10 [DDT programme ID 
279.B-5004(A), P.I. E. Dalla Bonta]. The NTT mounted the ESO Multi-Mode Instrument 
(EMMI) in red medium-dispersion spectroscopic (REMD) mode, using the grating No. 7 
with 600 grooves mm~^ in first order with a 1.0 arcsec x 5.5 arcmin slit. The detector was 
the No. 62 MIT/LL CCD with 2048 x 4096 pixels of 15 x 15 //m^. After a 2 x 2 on-chip pixel 
binning it yielded a wavelength coverage between about 4546 A and 6096 A with a reciprocal 
dispersion of 0.83 Apixel"^. The instrumental resolution was 2.5 A (FWHM) corresponding 
to ainst ~ 55 km s""*^ at 5900 A. The spatial scale was 0.332 arcsec pixel"^. We took three 
spectra of 2400 s along the major axis of the galaxy (PA= 54° ). The integration time was 
split into three exposures to deal with cosmic rays. Spectra of six giant stars with spectral 
type ranging from G7III to K5III were obtained, to use as templates. Arc lamp spectra 
were taken before and after every exposure to allow an accurate wavelength calibration. The 
average seeing FWHM during the observing runs was 1.0 arcsec in May 20 and 1.1 arcsec in 
July 10, as measured by the ESO Differential Image Motion Monitor. 

All spectra were bias subtracted, flatficld corrected, cleaned of cosmic rays, corrected 
for bad pixels and columns, and wavelength calibrated using standard IRAF routines. It 
was checked that the wavelength rebinning was done properly by measuring the difference 
between the measured and predicted wavelengths (Osterbrock et al. 1996) for the brightest 
night-sky emission lines in the observed spectral ranges.The spectra were co-added using the 
center of the stellar continuum as reference. The contribution of the sky was determined 
from the outermost ~ 30 arcsec and ~ 60 arcsec at the two edges of the resulting spectrum, 
respectively, and then subtracted. A one-dimensional sky-subtracted spectrum was obtained 
for each kinematical template star. Flux calibration was not performed. 

The stellar velocity dispersion was measured from the galaxy absorption features present 
in the wavelength range including the Mg line triplet (AA 5164, 5173, 5184 A) using the Fourier 
Correlation Quotient method (Bender 1990) as done in Corsini et al. (2007). IIR7429 (K3III) 
was adopted as the kinematical template. The stellar velocity dispersion was measured 
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within an aperture of l/Srg. It was then apphed the correction of Jorgensen et al. (1995) 
to derive the velocity dispersion within a circular aperture of radius l/8re. The correction 
was negligible, because the stellar velocity dispersion profile shows constant values, within 
the errors, toward the center. The stellar velocity dispersion of Abell 1836-BCG within a 
circular aperture of size l/8re (Tabled]) resulted to be cTc = 309 ± 11 km s~^. 

A consistent value (a^ = 308 ± 10 km s~^) was obtained by measuring the velocity 
dispersion with the Penalized Pixel-Fitting Method (Cappellari & Emsellem 2004) as done 
in Mendez-Abreu et al. (2008). The measurement of the velocity dispersion from spectra 
obtained in the Mg triplet region can be problematic for massive ellipticals, due to possible 
mismatch of the abundance ratios with those of the adopted stellar templates (e.g., Barth et 
al. 2003). To address this issue, the Mg region was masked out (Fig. 1201) . A lower value of 
= 288 ± 9 km s~^ was found after applying the aperture correction, and this was adopted 
for Abell 1836-BCG. 



6. Dynamical modeling 

Modeling of gas dynamics from HST data has steadily improved since the first successful 
applications of the method (Harms et al. 1994; Ferrarese et al. 1996; Barth et al. 2001; 
Sarzi et al. 2001; Marconi et al. 2003). For this study, we will use the procedure described 
in Coccato et al. (2006), where a detailed description of the method can be found. Briefly, 
a synthetic velocity field is generated assuming that the ionized gas is moving in circular 
orbits under the combined gravitational potential of stars and SBH. The gas is assumed to 
be confined in an infinitesimally thin disk centered at the photometric center of the galaxy. 
The model is projected onto the plane of the sky for a given inclination of the gaseous disk, 
and then degraded to simulate the actual setup of the spectroscopic observations. The latter 
step includes accounting for the width and location (namely position angle and offset with 
respect to the center) of each slit, instrumental PSF and charge bleeding between adjacent 
pixels. The free parameters of the model are the mass M, of the SBH, the mass-to-light 
ratio (M/L)^ of the stellar component (which we give in the J— band), and the inclination 
i of the gaseous disk; both {M/L)^, and i are assumed to be radially invariant. Although i 
can be estimated from the minor-to-major axial ratio of the disk, which is easily measured 
from the images, previous studies (Ferrarese et al. 1996, 1999; Sarzi et al. 2001; Shapiro et 
al. 2006) have shown that slight warps are common. When the ionized-gas emission arises 
predominantly from the innermost region of the dust disk, it is therefore best to treat i as 
a free parameter. The location of the slits, as well as the surface brightness distribution of 
the ionized gas (derived as described in § 13.51) are treated as input. M,, (M/L)*, and i are 
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determined by finding the model parameters which produce the best match to the observed 
velocity curve, obtained by minimizing = J^i"^ ~ v modY I ^"^ iv) where v ± 5{v) and Vmod 
are the observed and model velocity along the different slit positions, respectively. 

It has been noted in a number of studies of the kinematics of gas disks (Ferrarese et al. 
1996, 1998; van der Marel & van den Bosch 1998; Verdoes Kleijn et al. 2000, 2002, 2007) 
that the measured velocity dispersion of the emission lines greatly exceeds what is expected 
from instrumental broadening and radial variations in rotational velocity of the gas within 
a single slit element, and requires an intrinsic component of random motion. In our model, 
this is assumed to be described by a radial function of the form a{r) = + G\er'''l'^'^ . In 
principle, the exact characterization of a{r) affects the recovered velocity field, and therefore 
ctq, (Ji, and r„ should be treated as free parameters in the model, along with M, , i, and 
[M/L)^. In practice, however, previous studies (Verdoes Kleijn et al. 2000; Barth et al. 
2001; Coccato et al. 2006) have found that ctq, cti, and are fairly insensitive to the values 
of M. , i, and (M/L)^. 

We have run a set of preliminary models in order to find the shape of the intrinsic 
velocity dispersion profile that is required to match our data, within the typical range of 
our key model parameters. Specifically, we explored a reasonable grid of values for black 
hole mass, inclination, and stellar /—band mass-to-light ratio, finding for each combination 
of these parameters the values of (Tq, o"!, and that minimized = "^{v — VmodY / ^'^{v) + 
^((T — o"mod)^/^^(o") where a ± 5{(y) and CTmod are the observed and the model velocity 
dispersion along the major axis. In our final modeling procedure we have then fixed (Jq, 
(Ji, and Tcr to the values that were found for the best of such preliminary models (see next 
Sections) . 

In the following, the results for each galaxy will be discussed separately. 



6.1. Abell 1836-BCG 

We explored a three-dimensional grid of models with < M, < 3.6 x 10^° M© in 

2.4x10^ Mq steps, 0° < i < 90° in 1° steps, and < (M/L)^ < 33 (M/L)0 in 0.4 (M/L)© 
steps. Within this range, the intrinsic velocity dispersion of the gas best able to reproduce 
the obscrvables is a{r) = 12 + 204 e"''/'^'^ '^^'^ km s~^. The model adopts three parameters to 
fit 28 data points (the values of the rotational velocity), for a total of 25 degrees of freedom. 
The best model fitting the observed rotation curve (given the intrinsic velocity dispersion 
profile given above) requires M. = ?>Mtofo x 10^ Mq, i = 69.0°l}:^o, and {M/L)^= b.OtH 
(M/L)0(in the /—band), where the errors on M,, i, and (M/L)^, are quoted at the la 
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confidence level. The model, which has = 85.3 and a reduced Xr = 3.4, is compared to 
the observed spectrum and kinematics in Figure [21] and Figure [221 respectively. 

The small-scale asymmetries present in the model rotation curve (most evident at 
the central position) are a direct consequence of using in the models the intrinsic, PSF- 
deconvolved H«+[NII] intensity map (bottom panel of Figure [22]) . The observed rotation 
curves are well reproduced in general; the slight mismatch observed at r ~ 0".l for the slit 
located on the North- West side (-O'fl) of the nucleus (Figure [2^ might be due to the pres- 
ence of a slight warp in the outer part of the gas disk, which is not reproducible within the 
model's framework. 

The inclination angle derived from modeling the kinematic data {i = 69.0°li;9o), which 
extend only as far as r < 0".2, is significantly different from the inclination angle suggested 
by the axial ratio of the outer edge (r ^ 0".5) of the disk, i = 46° ± 5°. This suggests 
the presence of a warped structure, although such a tilting does not need to extend to the 
innermost regions where the kinematics was measured, and effectively our images can not 
constrain the gas geometry. 

Figures [231425] show the la, 2a and 3a confidence levels (following the Ax^ variations 
expected for two free parameters, i.e., 2.30, 6.17 and 11.8, Press et al. 1992) in the two- 
dimensional space of two of the fitted parameters, where the third parameter is held fixed 
at the best fitted value listed above. Finally, Figure [22] shows la, 2a, and 3a confidence 
levels individually on M,, i, and {M/L)^,, according to the Ax^ variations expected for one 
parameter (i.e., 1, 4, and 9, Press et al. 1992), marginalizing over all other parameters. 



6.2. Abell 3565-BCG 



The rotation curves along the three slit positions, for a total of 90 velocity data-points, 
were fitted for a grid of model parameters defined by < M, < 6.4 x 10^ Mq in 2.6x10^ 
M© steps, 0° < i < 90° in 1° steps, and < (M/L)^ < 12 (M/L)0 in 0.3 (M/L)0 steps. 
Within this range, the intrinsic velocity dispersion of the gas best able to reproduce the 
observables is a{r) = 122 + 98e~''/'^'°^ km s~^. The best model, which has = 1518 
and a reduced Xr = ^'''•5 (for 87 degrees of freedoms), requires M, = 1.34l'^Q;^g x 10^ Mq, 
i = 66.0°l||o and {M/L)^= 6.31^;^ (M/L)©, where all errors are given at the la confidence 
level. Such formal uncertainties were conservatively estimated after scaling up our errors on 
the velocity measurement until x^ = 1 (following Sarzi et al. 2001), acknowledging that our 
thin-disk model fails to reproduce some of the velocity structure observed in particular at 
the central and South-Eastern slit positions. Contrary to the case of Abell 1836-BCG, the 
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dynamical inclination angle of the disk is much closer to the morphological inclination angle, 
i = 71° ±1°, derived from the axial ratio of the outer edge of the disk (r ^ 1"A). Although our 
kinematic data extend almost to this radius (r < 0''8) in practice the luminosity weighting 
assigns more leverage to the data from the innermost region. The model is compared to the 
observed spectrum and kinematics in Figure [27] and Figure [281 respectively. 

Figures [291431] show the la, 2a, and 3a confidence levels in the two-dimensional space 
of two of the fitted parameters, where the third parameter is held fixed at the best fitted 
value listed above. Figure [52] shows la, 2a, and 3a confidence levels on M,, i, and (M/L)^ 
alone, marginalizing over all other parameters. It is worth noticing that mass-to-light ratios 
derived by the models refer to the nucleus of the galaxy, and might differ from the stellar 
population characterizing the galaxy on larger scale. 



6.3. Abell 2052-BCG 

In the case of Abell 2052-BCG, the asymmetric nature of the gas kinematics (§4.2) 
prevents the application of the method described above to pinpoint M,. An upper limit on 
the SBH mass can however still be obtained as described in Sarzi et al. (2002). In brief, an 
upper limit can be derived under the assumption that the central [N II] A6583 emission-line 
width, a, is produced by gas moving on circular orbits in a coplanar, randomly oriented disk 
in the Keplerian potential of the central black hole. If the line broadening is caused in part 
by non-gravitational forces, the derived upper-limit would well exceed the actual SBH mass. 

The model was implemented under the assumption that the intrinsic radial flux profile 
of the gas can be described by a Gaussian: 

Fir) =Fo + F, exp[-(r - ro)V24 J, (4) 

where the parameters Fq, Fi, tq, and afiux are derived by matching F(r), after convolution 
with the STIS PSF, to the observed emission flux profile along the slit. For a given black-hole 
mass, and under the previous assumptions regarding the gas geometry and kinematics, the 
model produces two-dimensional projected maps for the moments of the line-of-sight velocity 
distribution (LOSVD) at any position (x, y) on the sky: 

E^(x, y)= ! LOSVD(x, y, v:)v\^v, [k = 0, 1, 2). (5) 



Hv^ is then convolved with the STIS PSF; finally, the velocity dispersion, a, which cor- 
responds to (Sf2/Sf° — (Sw^/Sf °)^)^/^, is extracted within a 0.05" x 0.1" aperture for 
comparison to the observables. 



- 20 - 



There is no information on the orientation of the gaseous disk within the central aper- 
ture, nor is it possible to estimate it by analyzing the dust morphology, which is irregular 
and does not show a disk-like structure. In principle, a lower-limit on the SBH mass can 
be estimated assuming an edge configuration, whereas a model with a perfectly face-on 
disk would allow for an infinite SBH mass. Fortunately such extreme configurations are 
statistically rare, and considering that randomly oriented disks have uniformly distributed 
cosi we can estimate la upper and lower limits on M, by using models with i = 33° and 
i = 81°, corresponding to cosi = 0.84 and 0.16, respectively. For each i, the SBH mass was 
varied in 0.1x10^ M© steps between 6 x 10^ < M, < 11.0 x 10^ M© until the observed 
gas velocity dispersion was matched. The models were run for three values of mass-to-light 
ratio: (M/L)^=0, 5.0 (M/L)© (derived for Abell 1836-BCG), and 6.3 (as in the case of Abell 
3565-BCG) (M/L)q. The most conservative upper limit is obtained by neglecting altogether 
the stellar contribution to the gravitational potential, {M/ L)^=0, giving M,< 4.6 x 10^ Mq 
and M,< 1.8 x 10^ M0 for i = 33° and i = 81° respectively. Consistent with Sarzi et al. 
(2002), adopting (M/L)^=5.0 or 6.3 (M/L)©, has very little effect: M, decreases slightly to 
3.6 X 10^ M0 and 3.3 x 10^ respectively for i = 33° and 9.8 x 10^ and 7.7 x 10^ Mq 
respectively for i = 81°. Figure 1551 shows the dependence of M, on the assumed (M/L)^ for 
i = 33° and i = 81°. 



7. Discussion and conclusions 

We have presented a dynamical analysis aimed at constraining the masses of the SBHs 
in three brightest cluster galaxies: Abell 1836-BCG, at a distance of 147.2 Mpc, Abell 
2052-BCG, at a distance of 141.0 Mpc, and Abell 3565-BCG, at a distance of 50.8 Mpc. 
The models are based on data obtained with the Hubble Space Telescope. Broad-band 
WFPC2 and broad and narrow-band ACS images were used to constrain the luminosity 
profile and the distribution of the ionized gas, as traced by the Ha-|-[N II] emission, while 
STIS was employed to measure the rotation velocity and velocity dispersion profiles, from 
the [N II] A6583 emission, along three parallel slits positions, the first aligned along the 
photometric major axis, and the others adjacent to the first on either side of the nucleus. 

In the case of Abell 1836-BCG and Abell 3565-BCG, the regular morphology and 
kinematics observed for the ionized gas led to secure black hole mass detections of 
S.QI^q'^I X 10^ Mq and 1.34lQ;^g x 10^ Mq (where the uncertainties represent la errorCl), 



^Barth et al. 2001 show that the total error budget for gas-dynamical modelling is likely to be twice as 
large as such formal error estimates. 
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respectively. For Abell-2052-BCG, which displays irregular kinematics, an upper limit of 
M,< 4.60 X 10^ Mq was derived under the conservative assumption of a negligible stellar 
contribution to the gravitational potential, and an inclination angle for the gas disk of 33°. 
At face value, the black hole in Abell 1836-BCG is the most massive to have been dynamically 
measured to-date. 

In our modeling we have neglected the potential impact on our SBH mass determination 
of assigning a dynamical origin to the additional kinematic broadening that we have used 
in our models to match the observed profile for the gas velocity dispersion. Although it has 
been argued that this additional turbulence does not affect the bulk motion of the gas (e.g. 
van der Marel & van den Bosch 1998), it is important to notice that ignoring the possibility 
that the gas might be supported by dynamical pressure will lead to under-estimate the SBH 
mass. For instance, Barth et al. (2001) have shown that including a classical asymmetric 
drift correction led to an 12% increase of their best M, value. Unfortunately, we could not 
repeat this analysis in our models since the observed line-broadening exceeds by far the 
circular velocity, breaking down the epicyclic approximation on which the asymmetric drift 
correction is based. Nonetheless, our best SBH mass measurements are still likely to be 
underestimations of the real SBH mass. 

Figure [31 shows the location of Abell 1836-BCG, Abell 2052-BCG, and Abell 3565-BCG 
in the M,—ac (Ferrarese & Ford 2005; Tremaine et al. 2002) and near-infrared M,—Lbuige 
(Graham 2007; Marconi & Hunt 2003) planes. iT-band Two-Micron All-Sky Survey (2MASS) 
magnitudes were retrieved from the NASA/IPAC Infrared Science Archive, and corrected 
for Galactic extinction following Schlegel et al. (1998). Stellar velocity dispersions are from 
this paper (see g5]) for Abell 1836-BCG, from Smith et al. (2000) in the case of Abell 3565- 
BCG (we adopt his high S/N measurement, a = 335 ± 12 km s^^), and Tonry (1985) for 
Abell 2052-BCG. We note that published estimates of a for Abell 2052-BCG range between 
250 km s~^ and 370 km s~^; the high end of this spread is probably due to contamination 
from a companion galaxy located only 8" away in the North-East direction. The contribution 
from this galaxy was explicitly accounted for by Tonry (1985), whose value, a = 253 ± 12 
km s~^, we have adopted. Correcting the adopted values of a to a circular aperture of size 
l/8re (TableOQ), following Jorgensen et al. (1995), we find = 288±9 km s"\ = 233±11 
km s-\ and = 322 ± 12 km s"^ for Abell 1836-BCG, Abell 2052-BCG and Abell 3565- 
BCG, respectively. 

The SBH mass detected in Abell 3565-BCG, M,= l.SAtoH x 10^ M© , is consistent 
both with the M,— cXc (Ferrarese & Ford 2005) and the i^-band M,—Lhuige (Graham 2007) 
relations, which predict M,= 1.7^q;9 x 10^ M© (the error is computed adopting a 0.34 dex 
scatter in M.) and M,= l.l^oi ^ ^0 (adopting a 0.30 dex scatter), respectively. 
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For Abell 2052-BCG, although the conservative upper hmit on the black hole mass 
obtained for reasonable assumptions for (M/L)^ and i (§5.3), is consistent with both the 
M,—ac (Ferrarese & Ford 2005) and the K-hand M,—Lhuige (Graham 2007) relations, these 
predict significantly different masses, 3.5^^g x 10^ Mq (adopting a 0.34 dex scatter) and 
l-3^og X 10^ Mq (adopting a 0.30 dex scatter) respectively (but note the above caveat 
regarding the measurement of a for this galaxy). Although the sense of the discrepancy is 
consistent with the one argued for by Lauer et al. (2007), the lack of a precise determination 
of the SBH mass in this galaxy prevents us from establishing which, if either, of the two 
relations better predicts M, in this case. The only two BCGs in the hterature with a 
dynamically measured SBH mass are M87 (Macchetto et al. 1997) and NGG 1399 (Houghton 
et al. 2006). M87 appears to have a more massive black hole than expected based on its 
fT— band luminosity while it obeys the M,— (Jc relation. NGG 1399 is consistent with both 
the relations within their scatter. 

Abell 1836-BGG appears to be an outlier in all SBH scaling relations. The dynamically 
measured SBH mass reported in this paper, M, = 3.611q'5q x 10^ M0, is larger, at the 3a 
level, than the value of M,= 9. 614! ^ ^0 (adopting a 0.30 dex scatter), predicted by the 
i^T— band M, — Lhuige relation (Graham 2007). It is also larger, at the 3(T level, than the value 
of M,= 9.77^^.^3^ X 10^ Mq (adopting a 0.34 dex scatter), predicted by the M,— Uc relation 
of Ferrarese & Ford (2005). Finally, the SBH mass of Abell 1836-BGG is not consistent with 
the value of M, = 6.314^'^;^ x 10^ Mq predicted by the fundamental plane relation for SBHs 
by Hopkins et al. (2007), adopting an effective radius Tg = 13''1 = 9.3 kpc with a scatter of 
0.45 dex in M,. 

The existence of extremely massive SBHs as outliers to the empirical M,—ac relation 
would not be surprising. As shown by Lauer et al. (2007), the M,—Li,uige relation predicts 
a greater abundance of the most massive SBHs, given the luminosities of the BGGs and the 
lack of galaxies with velocity dispersions larger than ~ 350 km s~^. At the same time, the 
observed luminosity function of AGN implies the presence of SBHs as massive as ~ 5 x 10^ Mq 
within the distance of the clusters we have observed, and SBHs mergers can increase the 
highest SBH masses in these massive clusters up to ~ 10"*^° Mq (Yoo et al. 2007). A larger 
sample of SBH masses measured in massive clusters would be required to test models of the 
effect of mergers in increasing the largest SBH masses in the universe. 

In conclusion, both the M,—ac and M,—Lbuige relations appear at their high-M, end to 
be consistent with the SBH mass measured for one BCG, Abell 3565-BCG, but inconsistent 
with another one, Abell 1836-BCG. For the remaining target, Abell 2052-BCG, although the 
ionized-gas kinematics allowed us only to set an upper-limit on M, it seems unlikely that 
this galaxy could obey both relations simultaneously. The fact that Abell 1836-BGG is an 
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outliers in both the M,—ac and M,—Lbuige relations would appear to weaken the claim of 
Lauer et al (2007) that Lhuige is more reliable predictor of M, for BCGs. Overall, our results 
might indicate that the scatter of SBH scaling relations increases at the high end, although 
additional data are necessary to test this claim. 
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Table 1: Basic parameters of the sample galaxies 



Galaxy 


Other name 


Type 


1^25 


re 


Ills 


mx 


CZ 


D 


CTc 






[RC3] 


[arcmin] 


[arcsec] 


[mag] 


[mag] 


[kms-i] 


[Mpc] 


[kms-1] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Abell 1836-BCG 


PGC 49940 


SAO": 


1.4 


13.10'' 


14.56 ±0.14* 


9.99 ± 0.04 


11036 


147.2 


288 ± 9*^ 


Abell 2052-BCG 


UGC 9799 


E 


1.8 


55.98*' 


14.4 ± O.OS*" 


9.55 ± 0.06 


10575 


141.0 


233 ±11*' 


Abell 3565-BCG 


IC 4296 


E 


3.4 


41.4ia 


11.61 ± 0.05^^ 


7.50 ± 0.02 


3834 


50.8 


322 ± 12<= 



Note. — Col. (1): Galaxy name from Abell (1958) and Abell et al. (1989). Col. (2): Alternative 
identification. Col. (3): Morphological classification from RC3. Col. (4): Apparent major isophotal 
diameter measured at the surface brightness level \xb = 25 mag arcsec"^ from RC3. Col. (5): Effective 
radius (a: from RC3, b: from Hudson ct al. 2001). Col. (6): Total apparent B— band magnitude (a, c: from 
RC3, b: from Postman & Lauer 1995). Col. (7): Total apparent if— band magnitude from 21VIASS. Col 

(8) : Systemic velocities in the frame of the Cosmic Microwave Background, from Laine et al. (2003) Col. 

(9) : Distances from systemic velocity and Hq = 75 km s~ Mpc~^, except for Abell 3565-BCG (Jensen et al. 
2003). Col. (10): Stellar velocity dispersion (a: this paper; b: Smith et al. 2000; c: Tonry 1985) corrected 
to a circular aperture of size l/8re. 
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Table 2: Log of the ACS and WFPC2 observations. 



Galaxy 


Filter 


Pivot A 


Width 




Exp. Time 


Obs. Date 






A 


A 








(1) 


(2) 


(3) 


(4) 




(5) 


(6) 


ACS 


Abell 1836-BCG 


F435W 


4311 


310 




2 X 895.0 


2004 Feb 20 




FR656N 


6803 


136 




2 X 1382.5 


2004 Feb 20 


Abell 2052-BCG 


F435W 


4311 


310 




2 X 900.0 


2003 Sep 3 




FR656N 


6788 


136 




2 X 1380.0 


2003 Sep 3 


Abell 3565-BCG 


F435W 


4317 


293 




2 X 900.0 


2004 Aug 7 




FR656N 


6645 


133 




2 X 1388.0 


2004 Aug 7 


WFPG2 


Abell 1836-BCG 


F814W 


8012 


1539 




2 X 500.0 


2001 Apr 11 


Abell 2052-BCG 


F814W 


8012 


1539 




2 X 500.0 


2001 Apr 17 


Abell 3565-BCG 


F814W 


8012 


1539 


1 X 900.0 + 5 x 1100.0 + 4 X 1300.0 


1996 Jan 19-20 



Note. — Col. (1): Galaxy name. Col. (2): Filter name. Col. (3): Pivot wavelength (ACS: Sirianni et al. 
2005; WFPC2: Heyer et al. 2004). Col. (4): Filter width (ACS: Sirianni et al. 2005; WFPC2: Heyer et al. 
2004). Col. (5): Total exposure time. Col. (6): Observation date. 
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Table 3: Log of the STIS observations. 



Galaxy 


Offset 


PA 


Exp. Time 

N 

(4) 


Obs. Date 


(1) 


n 

(2) 


[1 
(3) 


(5) 


Abell 1836-BCG 


0.0 
-0.1 
+0.1 


-127.9 
-127.9 
-127.9 


2040 + 2800 
2 X 2800 
2 X 2800 


2004 Feb 21 
2004 Feb 21 
2004 Feb 21 


Abell 2052-BCG 


0.0 
-0.1 
+0.1 


51.0 
51.0 
51.0 


2030 + 2800 
2 X 2800 
2 X 2800 


2003 Sept 2 
2003 Sept 2 
2003 Sept 2 


Abell 3565-BCG 


0.0 
-0.2 
+0.2 


75.4 
75.4 
75.4 


2010 

2735 
2735 


2003 Aug 16 
2003 Aug 16 
2003 Aug 16 



Note. — Col. (1): Galaxy name. Col. (2): Measured offset of the slit center with respect to the galaxy 
nucleus. Col. (3): Position angle of the slit, measured from north to east. Col. (4): Total exposure time. 
Col. (5): Observation date. 
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Tabic 4. [N II] A6583 kinematics. 



r 


V 


u 


r 


V 


a 




r 


V 


c 


I J 


[Km s J 


[Km s J 


I J 


[kms-l] 


[kms" 




I J 


[km s ^] 


rVvn o — ll 

[Km s J 










(5) 


(6) 




\ ' ) 


\°) 












Abcll 1836-BCG 












+0"1 offset 






major axis 








-0".l offset 




-0.23 


11382± 25 


80 ± 28 


-0.23 


11500 ± 28 


165 ± 27 


-0.23 


11418 ± 24 


133± 24 


-0.18 


11398± 51 


256 ± 52 


-0.18 


11422 ± 31 


289 ± 31 


-0.18 


11361 ± 17 


137± 18 


-0.13 


11402± 53 


369 ± 51 


-0.13 


11455 ± 18 


266 ± 22 


-0.13 


11341 ± 18 


200± 18 


-0.08 


11288± 26 


306 ± 27 


-0.08 


11413 ± 11 


294 ± 13 


-0.08 


11283 ± 18 


275± 19 


-0.03 


11307± 19 


286 ± 20 


-0.03 


11369 ± 11 


370 ± 13 


-0.03 


11237 ± 30 


396± 38 


0.02 


11199± 17 


228 ± 17 


0.02 


11118 ± 13 


362 ± 17 


0.02 


11094 ± 35 


260± 31 


0.07 


11202± 43 


346 ± 45 


0.07 


11001 ± 13 


208 ± 17 


0.07 


11082 ± 30 


301± 43 


0.12 


11094± 21 


132 ± 22 


0.12 


10993 ±11 


114 ± 12 


0.12 


11024 ± 60 


241± 64 


0.17 


11091± 28 


164 ± 29 


0.17 


10978 ± 20 


123 ± 21 








0.22 


11014± 14 


51 ± 19 


0.22 


10982 ± 21 


92 ± 23 
















Abell 2052-BCG 












+0'^1 offset 






major axis 








-0'(1 offset 




-0.25 


10189± 27 


103 ± 28 


-0.30 


10172± 56 


209± 


56 


-0.30 


10100 ± 17 


83± 18 


-0.20 


10220± 39 


229 ± 37 


-0.25 


10156± 34 


179± 


31 


-0.25 


10011 ± 33 


213± 30 


-0.15 


10247± 17 


158 ± 16 


-0.20 


10130± 17 


183± 


17 


-0.20 


9979 ± 13 


193± 13 


-0.10 


10318± 16 


169 ± 16 


-0.15 


10160± 15 


176± 


23 


-0.15 


9912 ± 11 


102± 15 


-0.05 


10330± 8 


132 ± 9 


-0.10 


10154± 10 


171± 


16 


-0.10 


9908 ± 7 


186± 9 


0.00 


10329± 11 


169 ± 12 


-0.05 


10220± 9 


150± 


13 


-0.05 


9930 ± 7 


199± 9 


0.0 


10348± 15 


211 ± 15 


0.00 


10436± 18 


294± 


20 


0.00 


10075 ± 12 


168± 16 


0.10 


10312± 23 


236 ± 22 


0.05 


10494± 23 


190± 


28 


0.05 


10186 ± 25 


215± 29 


0.15 


10311± 16 


176 ± 15 


0.10 


10364± 19 


187± 


21 


0.10 


10242 ± 26 


162± 35 


0.20 


10311± 22 


132 ± 22 


0.15 


10251± 23 


199± 


21 


0.15 


10236 ± 16 


147± 16 


0.25 


10313± 21 


97 ± 21 


0.20 


10251± 16 


133± 


16 


0.20 


10242 ± 21 


191± 21 


0.30 


10328± 29 


194 ± 29 


0.25 


10263± 12 


83 ± 


13 


0.25 


10296 ± 15 


73± 16 


0.36 


10268± 22 


97 ± 23 


0.30 


10301± 30 


138± 


30 








0.41 


10459± 35 


109 ± 36 


0.36 


10293± 88 


159± 


89 








0.46 


10312± 14 


38 ± 19 


0.46 


10293± 32 


113± 


33 
















Abcll 3565-BCG 












+0'^1 offset 






major axis 








-O'^l offset 




-0.76 


4018± 20 


86± 23 


-0.81 


4011 ± 13 


62± 16 


-0.71 


4015± 26 


135 ± 27 


-0.71 


4010± 11 


80± 12 


-0.71 


4045 ± 24 


133± 25 


-0.61 


4054± 18 


86 ± 21 


-0.66 


4020± 11 


97± 12 


-0.66 


4073 ± 16 


113± 17 


-0.56 


3999± 11 


69 ± 13 


-0.61 


4046± 12 


115± 13 


-0.61 


4057 ± 14 


88± 16 


-0.51 


4006± 16 


112 ± 17 


-0.56 


4032± 9 


107± 9 


-0.56 


4064 ± 13 


100± 14 


-0.46 


4017± 15 


113 ± 16 


-0.51 


4057± 10 


132± 11 


-0.51 


4102 ± 16 


119± 17 


-0.41 


3987± 15 


120 ± 16 


-0.46 


4077± 12 


154± 12 


-0.46 


4101 ± 15 


117± 16 


-0.36 


3976± 13 


110 ± 14 


-0.41 


4073± 11 


142± 12 


-0.41 


4085 ± 14 


111± 15 


-0.30 


3970± 12 


97 ± 13 


-0.36 


4054± 13 


177± 12 


-0.36 


4099 ± 13 


143± 14 


-0.25 


3954± 8 


82 ± 10 


-0.30 


4028± 12 


181± 12 


-0.30 


4116 ± 16 


193± 16 


-0.20 


3939± 8 


107 ± 9 


-0.25 


4032± 13 


213± 14 


-0.25 


4125 ± 12 


206± 12 


-0.15 


3898± 9 


147 ± 10 


-0.20 


4016± 10 


224± 10 


-0.20 


4077 ± 9 


204± 10 


-0.10 


3866± 6 


115 ± 7 
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Table 4 — Continued 



r 


V 


(T 


r 


V 




r 


V 




cr 


["] 


[km s"-'] 


[kms-1] 


["] 


[kms-l] 


[kms^^] 


["] 


[km 




[km s"-*-] 


(1) 


(2) 


(3) 


{^) 


(5) 


(6) 


(7) 


(8) 




(9) 


-0.15 


3970± 8 


227± 9 


-0.15 


4068 ± 19 


226± 14 


-0.05 


3833± 5 


112 ± 6 


-0.15 


3971± 8 


234± 8 


-0.10 


4116 ± 15 


251± 15 


0.00 


3824± 4 


109 ± 4 


-0.10 


3926± 5 


206± 5 


-0.05 


4084 ± 17 


399± 18 


0.05 


3815± 4 


97 ± 4 


-0.05 


3874± 4 


191± 4 


0.00 


3830 ± 13 


451± 16 


0.10 


3804± 4 


99 ± 5 


0.00 


3818± 3 


155± 3 


0.05 


3673 ± 6 


254± 8 


0.15 


3746± 


15 


195 ± 15 


0.05 


3791± 3 


138± 3 


0.10 


3607 ± 6 


203± 8 


0.20 


3725± 


12 


167 ± 12 


0.10 


3762it 5 


164± 5 


0.15 


3590 ± 8 


171± 10 


0.25 


3681± 


16 


176 ± 17 


0.15 


3729± 8 


201± 8 


0.20 


3572 ± 14 


149± 15 


0.30 


3685± 


16 


134 ± 16 


0.20 


3686± 9 


204± 9 


0.25 


3599 ± 8 


181± 8 


0.36 


3669± 


17 


155 ± 17 


0.25 


3655± 11 


165± 11 


0.30 


3551 ± 11 


166± 11 


0.41 


3629± 


22 


163 ± 23 


0.30 


3652± 10 


150± 10 


0..36 


3551 ± 14 


179± 14 


0.46 


3670± 


19 


168 ± 20 


0.36 


3648± 9 


130± 9 


0.41 


3583 ± 23 
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Note. — Col. (l)-(3): Galactocentric distance, line-of-sight heliocentric velocity (uncorrected for 
inclination), and linc-of-sight velocity dispersion (uncorrected for instrumental velocity dispersion) mea- 
sured from the [NII]A6583 Une in the -|-0'(1 offset spectrum. Col. (4)-(6); As in Col. (l)-(3) but for the 
major-axis spectrum. Col. (7)-(9): As in Col. (l)-(3) but for the — O'^l offset spectrum. 
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Fig. 1.— ACS/HRC images of Abell 1836-BCG. Tlie panels sliow the F435W full frame (left) and 
a zoom toward the center (center). The same central section is shown in the right panel for the 
FR656N image. The grey-scale used for the left panel (full frame) is kept the same in Figures [!]■ 
[3l while the gray-scale for the other panels changes from filter to filter in order to highlight the 
distribution of the nuclear dust. All images have been background subtracted. 
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Fig. 2.— ACS/HRC images of Abell 2052-BCG. See Figured] for details. 
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Fig. 3.— ACS/HRC images of Abell 3565-BCG. See Figure ffl for details. 
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Fig. 4.— WFPC2/F814W images of Abell 1836-BCG {left panels), Abell 2052-BCG {middle 
panels), and Abell 3565-BCG {right panels). The top panels show the full frame, while the bottom 
panels show a zoom toward the center. The grey-scale used for the top (full frame) images is 
kept the same for all the three galaxies, while the gray-scale for the central sections changes from 
galaxy to galaxy to highlight the distribution of the nuclear dust. All images have been background 
subtracted. 
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Fig. 5. — Isophotal parameters for Abell 1836-BCG plotted against the "geometric mean" ra- 
dius r^eo = ~ with a measured along the semimajor axis of the galaxy. The panels 
show the radial profiles of: surface brightness (in mag arcsec"^) in ACS/F435W (squares) and 
WFPC2/F814W (filled circles) bands, F435W-F814W (filled circles) color, ellipticity e, position 
angle 6 (in degrees, measured from north to east), and parameters A3, A4, Bs, and B4, measuring 
deviations of the isophotes from pure ellipses (see text for details). Open circles show the inner 
radial profile of fi measured on the WFPC2/F814W image corrected for dust obscuration (as de- 
scribed in the text) with the addition of the Galactic extinction correction by Schlegel et al. (1998). 




Fig. 6.— 



The same as for Figure El but for Abell 2052-BCG. 




Fig. 7.- 



The same as for Figure [5l but for Abell 3565-BCG. 



- 39 - 




Si 



Abell 1836-BCG 



J 1 



0.0 0.2 0.4 0.6 0.8 1.0 

A(V) [mag] 



Fig. 8. — Optical depth map of Abell 1836-BCG (/e/i panet), and continuum-subtracted emission- 
line images in the ACS/FR656N band pass before {central panel) and after deconvolution {right 
panel). The field of view is 6".9 x 6''9. Orientation and scale are given in the left panel and they 
are kept the same for all the panels. 
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Fig. 9.— As in Figure El but for Abell 2052-BCG. 
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Fig. 10.— As in Figure El but for Abell 3565-BCG. 
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Fig. 11. — Deprojection steps for the stellar mass profile of Abell 1836-BCG. Top panel: Multi- 
Gaussian fit to the PSF-deconvolved surface brightness profile derived from the extinction-corrected 
WFPC2/F814W image. Middle panel: Deprojected stellar luminosity density profile. Bottom panel: 
Circular velocity curve assuming {M/L)^= 1. 
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12.— As in Figure CH but for Abell 2052-BCG. 
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Fig. 13.— As in Figure CH but for Abell 3565-BCG. 




Fig. 14.— a) HST STIS/F28X50LP acquisition image of Abell 1836-BCG. The image has been 
rotated to the STIS instrumental frame. Orientation and scale are given. The white cross is the 
position of the nucleus from STIS acquisition procedure. The rectangles overplotted on the image 
show the actual locations of the slit during the spectroscopic observations, b) Portions of the bidi- 
mensional STIS spectra obtained in position 1 (major axis of the dust disk, central panel), 2 {—0".l 
offset, left panel) and 3 (+0"1 offset, right panel). The spectral region centered on the Ha emis- 
sion line is shown after wavelength calibration, flux calibration and geometrical rectification. The 
spatial axis is horizontal and ranges between — 0''5 and +0'!5, while the wavelength axis is vertical 
and ranges from 6795 to 6856 A. Individual emission lines ([N II] A6548, Ha, and [N II] A6583) are 
identified on the right side of the figure, c) \N II] A6583 kinematics from the spectra obtained in po- 
sition 1 (central panels), 2 (right panels) and 3 (left panels). For each slit position the line-of-sight 
velocity curve (top panel), the radial profile of the line-of-sight velocity dispersion (uncorrected 
for instrumental velocity dispersion, middle panel), and the radial profile of line flux in arbitrary 
units (bottom panel) are given. Heliocentric velocities and velocity dispersions are not corrected for 
inclination and instrumental velocity dispersion, respectively. 
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Fig. 15. — As in Figure [T4l but for Abell 2052-BCG. In b) the spatial axis is horizontal and ranges 
between — 0''5 and +0".5, while the wavelength axis ranges from 6746 to 6834 A. 
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Fig. 16. — As in Figure [T^ but for Abell 3565-BCG. In b) the spatial axis is horizontal and ranges 
between —1".0 and +1".0, while the wavelength axis ranges from 6608 to 6697 A. Positions 2 and 3 
correspond to a slit-position offset of —0".2 and +0"2, respectively. 
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Fig. 17. — Nuclear spectrum {upper panel) and spectrum extracted at r = +0''15 (lower panel), 
along the major axis of the dust disk of Abell 1836-BCG. Each spectrum is plotted with the 
Gaussian-model narrow (dotted lines) and broad (dashed lines) components, reconstructed blend, 
and residuals. 
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18. — As in Figure [171 but for Abell 2052-BCG. The outer spectrum was extracted at 
+0".W. 



- 49 - 



Abell 3565-BCG 




6600 6700 6800 6900 

A [A] 



19. — As in Figure [171 but for Abell 3565-BCG. The outer spectrum was extracted at 
+0'^25. 
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Fig. 20. — Best fit of the broadened spectrum of the template star (thick red fine) to 
the spectrum of Abell 1836-BCG (thin black line) used to measure the galaxy velocity 
dispersion. The hatched regions mark the spectral ranges excluded from the fit to minimize 
the mismatch of the abundance ratios. The difference between the galaxy spectrum and 
broadedned template spectrum is shown in the bottom panel. 
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Fig. 21. — Comparison between the bidimensional observed (upper panels) and synthetic spectra 
{lower panels) for the best-fitting model of Abell 1836-BCG. Spectra obtained in position 1 {central 
panels), 2 {right panels), and 3 {left panels) are shown. In the models we added a random noise to 
mimic the actual S/N of the observed spectra. 
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Fig. 22. — Observed [N II] A6583 kinematics (filled circles) along witli tlie best-fitting model (solid 
lines) for the SBH mass of Abell 1836-BCG. The observed and modeled velocity curve (top panels), 
velocity dispersion radial profile (central panels), and flux profile (bottom panels) are shown for 
the slit along the +0''l offset position (left panels), the major axis (central panels), and the —0".l 
offset position (right panels) of the gas disk. The dotted lines correspond to a model obtained with 
M, = Mq, i = 69.0°, and (M/L)^= 5.0 (M/L)0(in the /-band). 
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Fig. 23. — Locus of points of equal values around the minimum Xm value for the model of 
Abell 1836-BCG. The la, 2a, and 3a confidence levels expected for two free parameters are shown 
in the M,-i plane, with {M/L)i,= 5.0 (M/L)© corresponding to the value. The grid of models 
explored is shown by the dots. 
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Fig. 24. — As in Figure \23\ but for the la, 2a, and 3a confidence levels expected for two free 
parameters shown in the M,-{M/ L)^, plane with i = 69° corresponding to the Xm value. 
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Fig. 25. — As in Figure \23\ but for the la, 2a, and 3a confidence levels expected for two free 
parameters shown in the i-{M/ L)i, plane with M,= 3.61 • 10^ corresponding to the Xm value. 
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Fig. 26.— distribution for Abell 1836-BCG as a function of M, (left), inclination (center), and 
[M/L)^ [right). The dotted horizontal lines indicate the confidence levels on the best fitting values, 
marginalizing over all other parameters. 
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Fig. 27.— As in Figure EH but for Abell 3565-BCG. 




Fig. 28. — Observed [N II] A6583 kinematics {filled circles) along with the best-fitting model {solid 
lines) for the SBH mass of Abell 3565-BCG. The observed and modeled velocity curve {top panels), 
velocity dispersion radial profile {central panels), and flux profile {bottom panels) are shown for 
the slit along the +0''I offset position {left panels), the major axis {central panels), and the —0".l 
offset position {right panels) of the gas disk. The dotted lines correspond to a model obtained with 
M, = M0, i = 66.0°, and {M/L)^= 6.3 (M/L)0(in the /-band). 
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Fig. 29. — Locus of points of equal values around the minimum Xm value for the modeling of 
Abell 3565-BCG. The la, 2a, and 3a confidence levels expected for two free parameters are shown 
in the M,-i plane, with {M/L)^= 6.3 (M/L)0 corresponding to the value. 
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Fig. 30. — As in Figure \29\ but for the la, 2a, and 3a confidence levels expected for two free 
parameters shown in the M,-{M/ L)^, plane with i = 66° corresponding to the Xm value. 
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Fig. 31. — As in Figure \29\ but for the la, 2a, and 3a confidence levels expected for two free 
parameters shown in the i-{M/ L)i, plane with M,= 1.34 • 10^ Mq corresponding to the Xm value. 
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Fig. 32.— Same of Figure [Ml but for Abell 3565-BCG. 
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Fig. 33.— Upper limit of SBH mass of Abell 2052-BCG as function of the (M/L)^ assumed, 
for an inclination of the gas-disk i — 33° and i — 81°. 
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Fig. 34. — Left panel: Location of the SBHs masses of our BCG sample galaxies (diamonds) with 
respect to the M,— Uc relation by Ferrarese & Ford (2005). We plot the SBH masses based on 
resolved dynamical studies of ionized gas (open circles), water masers {open squares), and stars 
{filled circles) from Table 2 of Ferrarese & Ford (2005). The dashed line represents the Tremaine 
et al. (2002) relation. Right panel: Location of our BCG sample galaxies with respect to the 
near-infrared (if-band) M^—Lhuigg relation by Graham (2007). Data are taken from Table 2 of 
Graham (2007) and symbols are as in the left panel. Masses and magnitudes of NGC 5252, NGC 
6251 and Cygnus A were adjusted to a distance obtained with Hq = 75 km s~^ Mpc~^. The dashed 
line represents the Marconi & Hunt (2003) relation. In both panels, we added the data relative to 
NGC 1399 (Houghton et al. 2006) and the dotted lines represent the la scatter in M,. 



